Volume Dependence of Bound States with Angular Momentum 



Sebastian Konig,^ Dean Lee,^'^ and H.-W. Hammer^ 

^ Helmholtz-Institut fiir Strahlen- und Kernphysik (Theorie) 
and Bethe Center for Theoretical Physics, Universitdt Bonn, 53115 Bonn, Germany 

^Department of Physics, North Carolina State University, Raleigh, NC 27695, USA 

(Dated: January 15, 2013) 

We derive general results for the mass shift of bound states with angular momentum > 1 in 
a finite periodic volume. Our results have direct applications to lattice simulations of hadronic 
molecules as well as atomic nuclei. While the binding of S-wave bound states increases at finite vol- 
ume, we show that the binding of P-wave bound states decreases. The mass shift for D-wave bound 
states as well as higher partial waves depends on the representation of the cubic rotation group. 
Nevertheless, the multiplet-averaged mass shift for any angular momentum £ can be expressed in a 
simple form, and the sign of the shift alternates for even and odd £. We verify our analytical re- 
sults with explicit numerical calculations. We also show numerically that similar volume corrections 
appear in three-body bound states. 

PACS numbers: 12.38. Gc, 03.65. Ge, 21.10.Dr 



Introduction. With recent advances in computational 
power and algorithms, the physics of numerous quantum 
few- and many-body systems can now be investigated 
from first principles. Lattice simulations are an impor- 
tant tool for such calculations [IHS] • The system is solved 
numerically using a discretized spacetime in a finite vol- 
ume. In practice, the finite volume is usually a cubic box 
with periodic boundaries. The box modifies the bound 
state wave functions and leads to shifts in the binding en- 
ergies. This shift needs to be subtracted from the calcu- 
lated energies for comparison to experiment. In the case 
of S-wave bound states, Liischer has derived a formula for 
the finite volume mass shift of two-body states [4 . See 
Ref. [5] for a recent application of this method in lattice 
QCD to extract the mass of the proposed H-dibaryon. 

But there are also many bound states with nonzero or- 
bital angular momentum. In nuclear physics, some par- 
ticularly interesting examples occur in halo nuclei [BHS]. 
These nuclei show a pronounced cluster structure. One- 
neutron halo nuclei can be regarded as a tightly-bound 
core with an extra neutron. In such cases the separation 
energy for the neutron is much smaller than the binding 
energy of the core as well as the energy required for core 
excitation. Thus, the volume dependence of energy lev- 
els obtained in ab initio lattice calculations of such halo 
systems would behave as a two-body system. 

A well known example of a P-wave halo state is the 
= l/2~ excited state in -'^^Be. The electromagnetic 
properties of the low-lying states in ^^Be can be well de- 
scribed in a two-body halo picture [H E]. If Coulomb 
interactions are included, proton halos like ^B also be- 
come accessible. In atomic physics, several experiments 
have investigated strongly-interacting P-wave Feshbach 
resonances in ^Li and [TU]-[T2] which can be tuned 
to produce bound P-wave dimers. There is interest in 
P-wave molecules in hadronic physics [13] as well as lat- 
tice investigations of the excited nucleon spectrum in a 



number of different spin channels [14) . Some of these 
states have been conjectured to have a molecular baryon- 
meson structure [15] . An extension of Liischer's formula 
to higher partial waves would provide a tool to discern 
molecular structures in hadronic states as well as halo 
structures in nuclei from the finite- volume dependence of 
lattice calculations for such systems. 

In this letter, we derive general formulas for the finite- 
volume mass shift for bound states with nonzero orbital 
angular momentum £. We also obtain a simple expres- 
sion for the multiplet-averaged mass shift for angular mo- 
mentum £. We verify our analytical results with numeri- 
cal calculations using an attractive short-range potential. 
We note recent studies on the related topics of extracting 
resonance properties and scattering phase shifts in higher 
partial waves from finite- volume energy levels [16|, I17[ . 
Although our analytic derivation can be applied rigor- 
ously only to two-body systems, we show numerically 
that quantitatively similar results also appear in three- 
body systems. 

Mass shift formula. In order to derive a general mass 
shift formula we consider a bound state solution lips) of 
the Schrodinger equation 

H\^b) = -Eb\^b) , H = -^Ar + Vir) (1) 

with angular quantum numbers {£, m) in a finite box of 
size with periodic boundary conditions. Following 
Liischer's derivation in [1], the energy shift compared to 
the infinite volume solution, 

AniB = Eb{(x) - Eb{L) , (2) 

can be written as 
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where n is an integer vector and k = ^J2iiEb is the bind- 
ing momentum. We assume that the potential has a finite 
range i? < L. Eq. (|3| arises from the overlap between 
copies of the system introduced by the periodic bound- 
ary conditions. For r > R, the wave function has the 
asymptotic form 
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where fi^ is a Riccati-Hankel function. We will use the 
relation 
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which follows from Lemma B.l in [18] and a derivative 
formula for spherical Hankcl functions [T^. The func- 
tions i?™ are the sohd harmonics defined via i?™(r) = 

For S- waves, Eq. (§ is a trivial identity. Inserting the 
Schrodinger equation to rewrite V{r) in (|3|, and using 
the fact that exp(— Kr)/(47rr) is a Green's function for 
the operator [A^ — k^] , we recover Liischer's result for 
S-wave bound states. In our notation this reads 
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We now generalize the mass shift formula to higher or- 
bital angular momentum. In the following we explicitly 
show the derivation for P-waves. We insert the asymp- 
totic expression Q with £ — 1 into Eq. ([S]) and use ([S]) to 
rewrite the Riccati-Hankel function. For m = we find 



Am 



(1,0) 



— T.Q-z'^Bir-nL) 

|n| = l 



r=0 



-0(e 



-V2kL\ 



(7) 



after integrating by parts. For m = ±1, the result is 
similar and involves derivatives with respect to x and y. 
Evaluating the sums then yields 
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Compared to the S-wave case, the P-wave mass shift is 
opposite in sign but equal in magnitude. Qualitatively, 
this means that S-wave bound states are more deeply 
bound when put in a finite volume while P-wave states 
are less bound. 

For higher partial waves we can proceed in exactly the 
same manner. The results, however, are more compli- 
cated and the shift for D-waves and higher partial waves 
depends on the quantum number m. We note that due to 
the periodic boundaries, the rotational symmetry group 



is reduced to a cubic subgroup. As a consequence, angu- 
lar momentum multiplets are split into irreducible repre- 
sentations of this subgroup (see, for example, Ref. [I6]). 
A similar splitting also arises in lattice calculations due 
to discretization artifacts. 

The mass shift for general £ can be expressed as 
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where the coefficients a are given in Tab. |l] for 

£ — 0,1,2. The irreducible representation of the cubic 
group is denoted by F in Tab. |Tj A detailed derivation 
of the general mass shift formula will be provided in a 
forthcoming publication [20 . The expressions for the fi- 
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TABLE L Coefficient a{x) in the expression for the finite vof- 
ume mass shifts for £ = 0, 1, 2. F indicates the corresponding 
representation of the cubic group. 

nite voiume mass shift become simpier when we sum over 
afl m for a given £. Using the trace formula for spherical 
harmonics, it can be shown i20i that 
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Dividing by 2i? -f 1, we obtain the average mass shift for 
states with angular momentum £. Apart from the overall 
sign, this average shift is independent of £. This follows 
from the fact that Yp{9, (j))Y(^*{6, (j)) averaged over m = 
—£, . . . ,£is equal to l/(47r) for all 0, (j), and £. For the case 
£ = 2 (cf. Tab. |l]), Eq. ^ can be verified explicitly by 
averaging over the three-dimensional representation 
and the two-dimensional representation . The mass 
shifts for the S- and P-wave states are especially simple 
because these multiplets are not split apart into more 
than one cubic representation. 

The sign of the finite volume mass shift can be ex- 
plained in terms of the parity of the wave function. At 
infinite volume the tail of each bound state wave function 
must vanish at infinity. At finite volume, however, the 
bound state wave functions with even parity along a given 
axis can remain nonzero everywhere. Only the derivative 
needs to vanish at the boundary, and the kinetic energy 
is lowered by broadening of the wave function profile. On 
the other hand, a wave function with odd parity along a 
given axis must change sign across the boundary. In this 
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case the wave function profile is compressed and the ki- 
netic energy is increased. We have illustrated both cases 
for a one-dimensional square well potential in Fig. [l] 





r\ 

/ 

1 — 1 — 1 




"^even 


^\ / 





X 



FIG. 1: (Color online) Wave functions with even (bottom) 
and odd parity (top) for a one-dimensional square well poten- 
tial in a box with periodic boundary conditions. The dashed 
lines give the infinite volume solutions for comparison. 

Comparison with numerical results. We test our pre- 
dictions for the S-wave and P-wave mass shifts using nu- 
merical lattice calculations. In Fig. [2j we show the mass 
shifts obtained from numerically solving the Schrodinger 
equation for a lattice Gaussian potential 

Vir) = -Vq exp (-rV(2i?2)) (11) 

with i? = 1, Vq = 6, and fi = 1. All quantities are in 
lattice units. This potential does not have a finite range 
in a strict mathematical sense, but the range corrections 
can be entirely neglected. In order to compare the depen- 
dence on the box size L with the predicted behavior, we 
have plotted log(L- lAmsl) against L (for S-waves Atob 
is negative). The expected linear dependence is clearly 
visible. 

For comparison we have calculated mass shifts using 
three different methods. The crosses show the direct 
difference, Eq. ([2|, where we have used L^o — 40 to 
approximate the infinite volume solutions. The boxes 
were obtained from the overlap formula ([s]). The circles 
were calculated using discretized versions of ^ and 
which we obtained by replacing exp(— Kr) /r with the lat- 
tice Green's function 

&(n) = G(„,-|) = iji;^5l^. (12) 

where (5^(q) = 2 J2i=i 2 3(1 ~ 1i) ■ "^^^^ Green's func- 
tion is also used to calculate the asymptotic normaliza- 
tion 7 from the lattice data. This incorporates the correct 
dispersion relation for our lattice model. Both the over- 
lap and Green's function results were calculated using 
lattice wave functions from the Lqo = 40 calculation. 

All three results agree well for both the S-wave and 
P-wave. For small L there are visible deviations which 



can be attributed to the ©(e"^"^) -corrections as well 
as potential range effects. The inset in Fig. [2] shows this 
more clearly. There we plot the relative differences be- 
tween the (logarithmic) direct results versus the overlap 
and Green's function data. 
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FIG. 2: (Color online) S- and P-wave mass shifts log(L • 
|Ams|) as functions of the box size L (in lattice units). We 
show the results obtained from the direct difference Eq. H 
(crosses), evaluation of the overlap integral Eq. (|3| (squares), 
and discretized versions of Eqs. (|6|, ([8| (circles). The dashed 
lines are linear fits to the squares. In the inset, we show the 
relative difference between the direct results and the overlap 
(squares) and Green's function (circles) data. 

When we perform a linear fit to the overlap integral 
data (dashed lines in Fig.[2| we obtain n — 2. 198 ±0.005, 
I7I — 11.5 ± 0.2 for the S-wave results; and k — 1.501 ± 
0.004, I7I = 7.0 ± 0.1 for the P-wave results. The values 
for the asymptotic normalization are in good agreement 
with the results I7I 11.5 (S-wave) and I7I 7.2 (P- 
wave) that are obtained directly from the Loo — 40 data. 
Inserting the corresponding energy eigenvalues into the 
lattice dispersion relation 

-hEb = (1 - cos(-iK)) , (13) 

we find K ^ 2.211 (S-wave) and k ^ 1.501 (P-wave), 
again in quite good agreement with the fit results. The 
remaining small discrepancies can be attributed to the 
mixing with higher partial waves induced by the lattice 
discretization and the fact that we have not performed a 
continuum extrapolation to vanishing lattice spacing. 

Summary and outlook. In this letter, we have derived 
an explicit formula for the mass shift of P- and D-wave 
bound states in a finite volume. We have compared our 
results with numerical calculations of the finite-volume 
dependence for a lattice Gaussian potential and found 
good agreement with predictions. For £ > 2 the mass 
shift depends on the angular momentum projection m 
due to different representations of the cubic group. The 
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average mass shift in a multiplet with arbitrary angular 
momentum (. can be expressed in a simple way, and apart 
from the alternating sign it is independent of £. Applica- 
tions to nuclear halo systems such as ^^Be and molecular 
states in atomic and hadronic physics appear promising. 
Our study provides a general framework for future lat- 
tice studies of molecular states with angular momentum 
in systems with short-range interactions. 

Finite-volume dependence can be used to probe the 
structure of a number of nuclei with conjectured alpha- 
cluster substructures piH^ . Recently, there have been 
ah initio lattice calculations of the low-lying states of ^^C 
using effective field theory [24]. In particular, the en- 
ergy for the spin-2 state of ^^C was calculated and found 
in agreement with the observed value of —87.72 MeV, 
just a few MeV below the triple-alpha threshold. It is 
not known how angular momentum is distributed in this 
state, and the study of finite volume effects may help to 
resolve this question. 

Our results apply rigorously only to two-body systems. 
However, there is empirical evidence that Eq. ([9| also 
gives the asymptotic L-dependence for three- and higher- 
body bound states at finite volume [53 dS]. In these 
cases the values of k and 7 are fitted empirically. We 
can show this explicitly using an extension of our Gaus- 
sian lattice model to three particle species. We take the 
particle masses to be equal, mi — 1112 = — 2, and 



consider Gaussian two-body potentials of the form (11) 
with the same range, R — 2, but different interaction 
strengths Fq^^ = 2.5, V^^ = 3.0, and V^^ = 3.5 between 
the particles 12, 23, and 31, respectively. 

In Fig. |3j we show the mass shifts for the lowest- 
lying trimer states with quantum numbers = 0+ and 
= 1^. All dimer states have less than half the binding 
energy of these trimers, which indicates that there is no 
underlying two-body molecular structure. As before we 










-5 






< 


-10 




-15 








-20 




-25 



'B, 



o .F - 



I I I I 

9 10 11 12 13 



L 



FIG. 3: (Color online) Trimer mass shifts log(— L- |Ams|) as 
functions of the box size L (in lattice units). We show lattice 
results for the two lowest-lying trimer states with j'^ — 0^ 
(squares) and = 1~ (circles). 

plot log(L- 1 Atob I) against L. Analogous to the two-body 
case, we find that Atub for — 0+ is negative, while 
Atub for = 1~ is positive. We also find the same 
linear dependence for log(L • lAmsl) at large L. Quite 



interesting are the subleading corrections which are es- 
pecially strong for = 1~ at smaller volumes. This 
may be due to competing finite volume corrections from 
negative S-wave and positive P-wave terms. These re- 
sults point to a possible new application of finite volume 
corrections to probe the radial distribution of angular 
momentum in complicated bound state systems. 
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